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Thermodynamic of universe with a varying dark energy component 
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We consider a FRW universe filled by a dark energy candidate together with other possible sources 
which may include the baryonic and non-baryonic matters. Thereinafter, we consider a situation 
in which the cosmos sectors do not interact with each other. By applying the unified first law 
of thermodynamics on the apparent horizon of the FRW universe, we show that the dark energy 
candidate may modify the apparent horizon entropy and thus the Bekenstein limit. Moreover, we 
generalize our study to the models in which the cosmos sectors have a mutual interaction. Our final 
result indicates that the mutual interaction between the cosmos sectors may add an additional term 
to the apparent horizon entropy leading to modify the Bekenstein limit. Relationships with previous 
works have been addressed throughout the paper. Finally, we investigate the validity of the second 
law of thermodynamics and its generalized form in the interacting and non-interacting cosmoses. 


I. INTRODUCTION 

Since the expanding universe is homogeneous and isotropic on scales larger than about 100-Mpc, it can be modeled 
by the so-called FRW metric [I| 


ds^ = dt^ — (t) 


dr^ 


1 — kr"^ 


■ r^dn^ 


( 1 ) 


where fc = 0, ±1 is the curvature constant corresponding to a flat, closed and open universe, respectively. Additionally, 

2 

a{t) is the scale factor written as a{t) = aotfor ui > —1 and a(t) = aoexpHt when uj = —1, whiles w = ^ is the 

state parameter of prefect dominated fluid. In addition, H = - is the Hubble parameter [ij. Moreover, for Phantom 
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regimes (w < —1) the scale factor is written as a{t) = ao{tbr — t) 3(i+“), where t^r is the big rip singularity time, 
everything will be decomposed to its fundamental constituents at that time, [^. Additionally, it is shown that one 
can use the conformal form of this metric to describe the inhomogeneity of the cosmos in scales smaller than 100-Mpc 
0. In the standard cosmology a primary inflationary expansion era is used to get a suitable theoretical description 
for horizon problem which emerges in the study of Cosmic Microwave Background (CMB) [ij. Observational data 
^nals us a universe with a > 0 and d > 0 which means that we need to modify the gravitational theory 

or considering an unknown source, named dark energy (DE), for describing this phase of expansion (lol - [l^ . A 
simple model used to explain DE considers an unknown fluid with constant density, pressure and wd = — 1 called the 
cosmological constant (CC), and leads to an exponential expansion (a(t) = aoexpHt) [l|. We should note that the 
current expanding phase of the universe is in full agreement with both of the thermodynamics equilibrium conditions 
and the rise of complexity content of the universe meaning that the universe may maintain its current expanding 
phase [l^. More studies on the thermodynamics of DE and the final state of universe can be found in refs. [T^ - [l^ . 
Bearing the primary inflationary era together with the CC model of DE in mind, two difficulties including the fine 
tuning and coincidence problems are inevitable [l|. It is also useful to mention here that since the CC model has a 
satisfactory match to the observational data, it formes a basis for the standard cosmology [l|. 

Observational data support a DE candidate with varying energy density p7l - l^ . Indeed, there are various attempts 
to model the source of the primary and current accelerating eras by introducing a varying model for the DE candidate 
[9l-[l^ [23l - l^ . Recently, Lima and co-workers proposed that a universe filled by a dynamical vacuum en erg y density 
can avoid the big bang as well as the big crunch singularities, the fine tuning and coincidence problems [8 (t| . Indeed, 
since the vacuum density is decreased as a function of the Hubble parameter, the Lima’s model has enough potential 
for solving the fine tuning and coincidence problems [Mill- Additionally, because in their model, the cosmos began 
to expand from a primary unstable de-Sitter spacetime, and finally reaches to another eternal de-Sitter spacetime, the 
horizon problem as well as the big crunch problem are naturally solved [s^. Moreover, It is shown that the ultimate 
de-Sitter spacetime is in accordance with the thermodynamic equilibrium conditions, and therefore the cosmos may 
serve its final stage 38|. In this model, the state parameter of the vacuum energy satisfies the wd = — 1 condition. 
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and thus the vacuum energy decays into the other fields confined to the apparent horizon of the FRW universe 
[ 3 ^ [ 3 ^ . It is worthwhile to mention here that the decay of vacuum into the other fields is due to a mutual interaction 
between the cosmos sectors leading to leave thermal fluctuations into the cosmos in this model [s^. It is a good 
feature, because observations allows a mutual interaction between the cosmos sectors [dol - l^ . Thermodynamics of 
such possible mutual interactions are also studied in various theories of gravity by considering various models for DE 
|44l - l4^ . In fact, the relation between such possible interactions, coincidence and fine tuning problems and thermal 
fluctuations attracts more investigators to itself [H, [H, |4^ . 

Similarity between the Black Holes laws and those of thermodynamics motivates us to define a temperature as 



( 2 ) 


where k is the surface gravity of Black Hole [s^. On one hand, for some spacetimes, such as the de-Sitter spacetime, 
surface gravity and thus the corresponding temperature are negative (50j |. and therefore we need to define T = ^ in 
order to get the positive values for temperature Whiles, on the other hand, one can get the Einstein equations 

on the event horizon of Black Holes (as a causal boundary) by applying the first law of thermodynamics on the event 
horizon and considering Eq. (I5|) as a suitable definition for temperature Indeed, it seems that this similarity 

is much more than a mere resemblance (53l - [6^ . 

The apparent horizon of the ERW universe, as the marginally trapped surface, is evaluated by 


daCd°‘C = 0^rH, 

where C = CL{t)r, and can be considered as the causal boundary Therefore, One gets l6i,|2| 


rA 


1 





( 3 ) 

( 4 ) 


Moreover, the surface gravity associated with the apparent horizon of the ERW universe can be evaluated by using 

^=-^da{V^h-’’dbC). (5) 

where hab = diag(—I, a(t)^) Since the WMAP data indicates a flat universe, from now we set k = 0 

Thus, simple calculations lead to 


and therefore 


k = -H{1 + 


H , 

2iJ2' 


( 6 ) 



H H , 
2^' 


( 7 ) 


where we have used Eq. ([^l to obtain this equation . It is useful to note here that for the FRW universe 

supported by a fluid with p = —y = constant {ui = —1), H = 0 and therefore this equation covers the result of 
de-Sitter spacetime (T = — ^) [50, [MI- 

In cosmological setups, some authors use various definition of temperature and get the corresponding Einstein equa¬ 
tions (Eriedmann equations) on the apparent horizon [SlL l52l . l68l . I7ll474l| . In order to avoid the negative temperature, 
authors in [s^, have defined T’ = ^ and used the first law of thermodynamics (in the TdSA = —dQ form) to 

get the Eriedmann equations. In their approach Sa = (the Bekenstein limit) and Q are the horizon entropy and 
energy crossed the apparent horizon, respectively (5l|. Indeed, authors argued that the extra minus sign in the first 
law of thermodynamics is the result of universe expansion leading to decrease the energy of confined fluid together 
with increase the size of the universe and thus Sa- Therefore, by using original definition of temperature ([2]) and 
thus o, called the Hayward-Kodama temperature [HI, [H, [ll, Q , together with the TdSA = dQ form of the first 
law of thermodynamics we can cover the Eriedmann equation. Moreover, it seems that W = ^ habT°-’’, where 
is the energy momentum tensor of fluid which s prea ds over the cosmos, plays the role of pressure in the dynamics 
spacetimes and thus the ERW universe [HI 1^. l7ll - [^ . Following this argument, authors in [HI, [S [HHzll have 
used (m and the work density definition (W) to get the Friedmann equations by applying the first law of thermo¬ 
dynamics {TdSA = dQ = dE — WdV) on the apparent horizon of the FRW universe in various theories of gravity, 
whiles E is the energy confined to the apparent horizon. It is also shown that Loop Quantum Gravity corrects the 
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horizon entropy which leads to modify the Friedmann equations on the apparent horizon if one considers Q together 
with TdSj\ = dE — WdV [H, [t^. The entropy of a self-gravitating system depends on the gravitational theory used 
to describe the gravity field. Accordingly, it seems that the self-gravitating systems satisfy the Bekenstein limit of 
entropy in the Einstein general relativity framework. But, since the origin of DE is unknown, it may have either a 
geometrical or physical origin, one can expect that the DE candidate may affect the horizon entropy. By the same 
token, it is shown that the ghost dark energy and its generalization, as the dynamics candidates for DE, may also 
add an additional term to the entropy of various horizons leading to modify the Bekenstein limit 1?^ . Therefore, 
it seems that the dynamics model of DE may lead to modify the horizon entropy and thus, the Bekenstein limit. 
Recently, it is also shown that a mutual interaction between the cosmos sectors may change the horizon entropy . 
The second law of thermodynamics states that the horizon entropy may meet the > 0 condition [^. Nowadays, 
thanks to the Bekenstein works [13, HH , it is believed that the rate of the total entropy of a gravitational system 
should be positive meaning that while Sin is the entropy of confined fluid. The latter is called the 

general second law of thermodynamics |5ll ISOl l8l| . Comprehensive reviews on the various temperature definitions in 
cosmological setups, their motivations to geth er with the validity of the first, second and generalized second laws of 
thermodynamics can be found in refs. [^1^1^. Now, one can ask how a DE candidate and its probable interaction 
with other parts of cosmos affect the horizon entropy, the second and generalized second laws of thermodynamics? 

In this paper, we point to the unified first law of thermodynamics and assume that it is available on the apparent 
horizon of the flat FRW universe, while T (the horizon entropy) corresponds to the Hayward-Kodama definition of 
temperature © on the apparent horizon of the FRW universe [6^, > and show that a DE candidate may lead 

to a new bound for the horizon entropy, whiles the cosmos sectors do not interact with each other. Additionally, we 
show that any mutual interaction between the cosmos sectors may also modify the horizon entropy. The relationships 
with similar works are also studied. Moreover, the results of considering the Cai-Kim temperature are also derived. 
Finally, the validity of the second law of thermodynamics and its generalization is also addressed. Since the physics 
behind the Lima’s model [s^ is completely different from the ordinary models, introducing for describing DE, we 
point to results of considering this model. 

the paper is organized as follows. In the next section, after a brief review on the previous related works, we apply 
the unified first law of thermodynamics on the apparent horizon of the flat ERW universe, and show that how a 
dynamic candidate for DE may change the horizon entropy, whiles the cosmos sectors do not interact with each 
other. Thereinafter, we generalize our study to the interacting case and get a modification for the horizon entropy 
due to the mutual interaction between the cosmos sectors. In section (III), we study the validity of the second law of 
thermodynamics and its generalization. Section (IV) is devoted to a summary and concluding remarks. Throughout 
this paper, we set G = c = h = 1 for simplicity. 

II. HORIZON ENTROPY AND THE UNIFIED FIRST LAW OF THERMODYNAMICS 

The unified first law of thermodynamics, which is available in some theories of gravity, is written as 

dE = A'li + WdV, (8) 

where W = ^habT^^ and A = §(1 - h^^da0b Ob=£. are the work density and the Misner-Sharp energy confined to 
the apparent horizon, respectively [b^, l76l - l78j |. In addition, A and are the area of horizon and the energy 

supply vector, respectively, and 


A'^ = Atjjadx°', 


while 


= T%C + WdaC, 


(9) 

( 10 ) 


is the projection of the total four-dimensional energy-momentum tensor in the normal direction of the two- 
dimensional sphere. Consider a perfect fluid source (T^ = diag{—pT,PT,PT,PT)) together with Friedmann equations, 
by simple calculations we get E = prV, 

dE - WdV = VdpT + + (11) 

and 




2 


2 


( 12 ) 
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where a is the scale factor. Using the energy-momentum conservation law = 0) 

Pt + 3i/(/0T + Pt) = 0, (13) 

and adr = dC — rda in rewriting Eq. (1121) to obtain 

A'if = VdpT + (14) 

where dV = AdQ and A = 47rC^. By comparing this equation with (fTTl) . we get 

A'H = dE- WdV, (15) 


which is the unified first law of thermodynamics. This result is independent of the number and nature of fluids which 
support the background spacetime. In addition, one may decompose T^i, into 


T 

fll/ 


rpDE 


+ T, 


flU 5 


(16) 


where and are the energy momentum tensors of DE and the material parts of cosmos (radiation, matter 
and etc.), respectively. In this situation, it is apparent that 4' = -|- dt™, W = + IE’” and E = E^^ + E™, 

where E^^ = and E™ = p^V. Therefore, by following the above argument, whenever = 0, 

we get 

AT™ = AT - AT^^ = dU™ - lE™dU. (17) 

l^^[(Uie^eat flow crossing the horizon) is determined by the pure matter energy-momentum tensor (T^) as [11,161- 

6Q = AT™. (18) 

Eor some gravitational theori^ one can use the Clausius relation together with Eq. 0 to get the horizon entropy 
{Sa) by using 

TdSA = SQ = AT™, (19) 


which leads to 


TdSA = at™ = at - AT^^ = dU™ - tE™dU, (20) 

where we have used dnl) to get the last equation. It is useful to note that for some theories such as the f{R) gravity, 
Eq. m and thus (uni) is not always available [^ . 

Recently, some authors considered the Hayward-Kodama definition of temperature 0, a universe filled by either 
a ghost dark energy or its generalized form together with a pressureless matter and use the TdSA = AT — AT^^ 
relation to get an expression for the entropy {Sa) [z 1|- Their results show that the entropy of the matter fields 
differs from the Bekenstein entropy due to the DE effects. They argued that their results are in agreement with the 
entropy of the apparent horizon in the DGP braneworld model which signals that this approach may be used to get 
the entropy of apparent horizon in other theories of gravity. This adaptation between these entropies signals that one 
may find a geometrical interpretation for the origin of the ghost dark energy model (as a DE candidate) by using the 
DGP braneworld model of gravity. Motivated by this work, Sheykhi extended their work to the apparent horizon of 
the ERW universe and used the TdSA = dE"^ — W^dV relation to get the same result as that of ref. [t^ for the 
entropy. Finally, he concludes that the obtained relation for the entropy (Sa) may be interpreted as the corrected 
relation for the apparent horizon entropy 0. It is useful to stress here that Eq. (HZl clarifies the reason of getting 
the same result for the horizon entropy by authors in Refs. [zilzl- Moreover, their results are available only when 
= 0 which means that the cosmos sectors do not interact with each other. Recently, by considering 
the ERW universe in which the cosmos sectors interact with each other, Mitra et al. use the TdSA = AT — AT^^ 
relation to get the trapping horizon entropy in the Einstein relativity framework. They argued that the obtained 
relation for the entropy differs from the Bekenstein entropy due to the mutual interaction between the cosmos sectors 
[^ . In continue, Mitra et al. extended their hypothesis to other different gravity theories [t^. Moreover, it is also 
shown that a gravitationally induced particle production process as the DE candidate may change the horizon entropy 
[ 3 ^. Bearing the Lovelock theory in mind, some authors used the TdSA = AT — AT^^ relation to get the entropy of 
the apparent horizon in cosmological setup . Another study including the loop quantum cosmology can be found 
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in ref . Here, by considering a varying DE candidate, we are going to find a general relation for the entropy of the 
apparent horizon in both of the interacting and non-interacting cosmoses and investigate the second and generalized 
second laws of thermodynamics in the Einstein relativity frame work where Eq. dm) and thus & are valid. 

For this propose, consider the flat FRW universe with Friedmann equation 

H^ = Y(p + Pd), ( 21 ) 

where pu is the density of dark energy component. In addition, p is the density of rest fluids in the cosmos which 
may include baryonic matters, dark matter and etc, leading to 

P = Pbm + Pdm + ■•■ (22) 

Therefore, p is nothing but p"* which is previously introduced. For the sake of simplicity, we omit the m label 
throughout the paper. From Eq. m and the Bianchi identity we get 

2HdH-^dpD = ^dp, (23) 


and 


p-f 37J(p-|-p)-|-p_D-|-3iJ(pD +Pd) = 0, (24) 

which is nothing but the energy-momentum conservation law, respectively. In this equation, pu and p denote the 
vacuum pressure and the pressure corresponding to the density p, respectively. Dot is also denoted as the time 
derivative. Consider a dark energy candidate with density profile 

3 a + + 37H’” 

PD =--, ( 25 ) 

which converges to CC whiles /3 = 7 = 0. Whiles n = 5 , it covers the ghost dark energy model and its generalization 
for CK = /3 = 0 and a = 0, respectively [23l - [^ . The 7 = 0 case has been extensively studied in the literatures [2^-[3lj. 
The results of considering either an arbitrary value for n or optional function of H for po = f{H) can be found in 

M- 

Moreover, the cosmological applications of considering model with a = 0, n = | and uju = — 1 has also been 
studied [11,133. More similar density profiles for the DE candidate with wd = — 1 can also be found in [Hi. Another 
attractive case proposed by Lima et al. is obtainable by imposing the n > 1 condition together with wd = — 1 to the 
density profile of the DE candidate, whenever n is also an integer number [H - It is useful to note that Ed = pnV 
and E = pV are the energy of dark energy component and the energy corresponding to the density p, respectively. 
Therefore, E is nothing but E"* mentioned previously and we omit the m label for the sake for simplicity. 


A. Non-Interacting Models 

At the first step we consider a universe in which the cosmos sectors do not interact with each other. Therefore, the 
energy-momentum conservation law implies (12411 


p + 3H{p + p) = 0, 

(26) 

and 


PD + 3H {pD +Pd) =0. 

(27) 

Substituting (1^ into (1^ to get 


^p = -a- - 777 ^". 

0 

(28) 

Bearing Eq. (|26ll in mind, by using Eq. (l23l) we reach at 


(2iJ(l — P) — 2n'yE['^'^~^)dH = —8TTH{p + p)dt. 

(29) 
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Using the Hayward-Kodama temperature relation (—T = ^(1 + 2 §^)) 01 obtain 

- T{2H{1 - /3) - 2n'jH^^-^)dH = -dH'^{p + p)dt - 2{p + p)dH. 

From Eq. (1^ . since E = pV and dV = —j^dH, we get dE = —ATvpE[~'^dH — AtvH~'^{p + p)dt leading to 


{p + p)dt = 


H^dE pdH 

IP' 


If we combine this equation with (1301) we obtain 


T{-2H{1 - P) + 2njH^'^-^)dH =- dE + 2{p - p)dH. 

7T 


It is easy to show that this equation can be rewritten as 

Ott 

r[(-^(l -/3)+ 2n77rij2"-5)d7J] = dE - WdV. 


(30) 


(31) 


(32) 


(33) 


In this equation W = is the work density required for applying a hypothetical displacement dr a to the apparent 
horizon [t^ . By comparing this result with Eq. (uni, one gets 


leading to 


O-TT 

dSA = {-j^{l-P) + 2n-fTTH^^-^)dH 


c _ ^ ,2-n 

SA--il-P) + ^^A , 


(34) 


(35) 


where A = ATrr\ = is the area of horizon. Therefore, ”^^2 ^ term besides the area term. In addition, 

since the entropy is not an absolute quantity, we have set the integral constant to zero. It is also apparent that, for 
n = 2, entropy is not well-defined. In order to eliminate this weakness, let us restart from Eq. ([Ml). by substituting 
n = 2 and taking integration from that, we get 


_ A yTTluTT 77rln^ 

OA — ^(1 - p)-^-1-^-1- c>o- 


(36) 


Finally, since entropy is not an absolute quantity, one can set So = and gets 

^ A, 77r In >1 , , 

Sa = -{1-P) + ^^. (37) 

Therefore, models with n = 2 induce a logarithmic correction to the horizon entropy. Logarithmic correction terms 
have been previously proposed by some authors which either consider the thermal equilibrium and quantum fluctu¬ 
ations in loop quantu m g ravity framework or the thermal fluctuations of system about its thermodynamic 

equilibrium state [93 1^- Indeed, logarithmic correction due to the thermal fluctuations are valid in all physical 
systems Let us study some choices with n = Bearing Eq. ([M]) in mind. For a constant vacuum energy density 
{pd = Oi), we face with the kCDM theory and we get Sa = ^ which is in agreement with previous studies [^ItiII^. 
Moreover, for a = 0, /3 = 0 and n = we have 


37 „ 


(38) 


which is the profile density of ghost dark energy model 


. In this limit, from Eq. (|M|) . we get 


<? =^- 
^ 4 3V^ 


7 .. 3 

' A^. 


(39) 


which is in agreement with the ghost dark energy modification to the entropy evaluated previously [ziliS- Here, we 
have used the Hayward-Kodama definition of temperature © together with the apparent horizon of the FRW universe 
to get this relation whiles, author in 771, has considered U = -^ to get (IMl) on the apparent horizon. Moreover, 
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authors in [t^ used trapping horizon and the temperature definition ^ to get this relation. Additionally, 

equation (ESI), for a = 0 and n = reduces to 


_ 3 ^ . 37 


(40) 


which is the profile density of generalized ghost dark energy model [25l . l82| . By considering this profile density we get 


as the modification of the generalized ghost dark energy model to the horizon entropy [t^. Although this result is 
previously obtained by authors in ref. (t^ . but our derivation is completely different from that of they. Here, we 
worked on the apparent horizon whiles they have considered the trapping horizon and found the similar results. In a 
more general case, for arbitrary functional form of pD^ by following the above recipe we get 


dSA “ 


(42) 


where prime denotes derivative with respect to H. Taking integral to obtain 


= 


A 

J 


3 




dpD + C, 


(43) 


where C is the integral constant. Therefore, a varying DE candidate imposes a correction term to the horizon entropy 
in accordance with the first law of thermodynamics and thus, the second term of the RHS of Eq. (l43l) . It is also useful 
to note that the result of considering CC (Sa = ^) is obtainable by substituting dpo = 0 in this equation [tiI. [ 73 |. 

Now, let us use the Cai-Kim temperature (T = j y ! to get the entropy of apparent horizon. In order to achieve 
this goal, we follow the approach of authors in ref. [^, where TdSA = —dQ and dV — 0. Using this argument and 
bearing Eqs. (HU) and (l20l) in mind to reach 


dSA = -^dp. 

Now, by substituting dp from Eq. (l23l) into this equation, one gets 


(44) 


which leads to 


dSA = 


27r 

IP 


dH 


Sa = 


A 


3 


Stt^ 

3IP 


dpD 


dpD, 


C, 


(45) 


(46) 


where C is the integration constant. Therefore, once again, we get a relation for the horizon entropy which is in full 
agreement with the previous result (H51) . obtained by considering the Hayward-Kodama temperature. 


B. Interacting Models 

When the cosmos sectors interact with each other, energy-momentum conservation law implies (j24[) 

p+ 3H{p + p) =-pD - iH{pD +Pd), (47) 

meaning that 

dp = —3H{p + p)dt — dpD — ^H{pd + pD)dt. (48) 

Therefore, by considering Eq. (l23l) and following the recipe which leads to Eq. (l4^ . we get 

dSA = -^dH - ^^{pD +PD)dt, 


(49) 





which yields 


S^ = ^-8n^ I P£^dt + C, 


(50) 


where C is again an integral constant. Therefore, the second term of RHS of this equation is nothing but the entropy 
correction due to the mutual interaction between the cosmos sectors. For interacting models in which the state 
parameter of the DE candidate meets the ujd = — 1 condition, and therefore pu + Pn = 0, this additional term is 
zero meaning that the horizon entropy in these models satisfies the Bekenstein limit [80|. For instance, in the model 
proposed by Lima et al. [3^ . in which vacuum decays into the other parts of cosmos and po + pd = 0, the horizon 
entropy of the flat FRW universe meets the Bekenstein limit [s^ . It is in agreement with the initial and final de-Sitter 
spacetimes of this model, since the horizon of de-Sitter spacetime meets the Sa = ^ condition [5^ [^, [^. Now, 
let us derive Eq. (I50|) by using the unified first law of thermodynamics. Bearing the definition of dt in mind, simple 
calculations lead to 


A^de ^ _ ^VipD+PD)H ^^ _ 


Mpd +Pd) 


[dC — CHdt], 


2 2 

where we have used the rda = C,Hdt relation to obtain this equation. It is a matter of calculation to show 




(51) 


(52) 


Eq. (fTg|l is valid, and thus, simple calculations lead to 


where we have used dV = —^Hdt to get this equation. Since we work in the Einstein general relativity framework, 


TdSA =A^- A^^^ = -^[1 + - ^(.PD+PD)dt), (53) 

where we have used the A^ = T{—^dH) relation, while T is the Hayward-Kodama temperature, in obtaining this 
relation [5l|, [bl, [t^, [t^, [ll| ■ It is apparent that this equation is nothing but (l4^ which leads to Eq. (15(1 . 

Our result is in agreement with the recent work by Mitra et al. [t^ . Whereas, we have started from the Friedmann 
equations and considered the apparent horizon as the causal bound, Mitra et al. used the trapping horizon and 
relation SQ^ = A'k — to obtain (l50l) . It is apparent that Eq. (l53l) clarifies that why both of us get the same 

results, while, our start points differ from each other. 

Finally, let us consider the Cai-Kim temperature to estimate the horizon entropy. In this situation, for an infinites¬ 
imal time dV = 0, and from Eqs. (lUl) and (ICT) we get 

TdSA = = -AT -h = -Vidp + dpo) - (54) 

where we have followed the approach of authors in ref. (5^ in order to define TdSA = —5Q. Now, bearing Eq. dan 
in mind, since T = ^, simple calculations lead to 

Sa = ^-8n^ J ^^^^dt + C, (55) 


where C is an integration constant. Therefore, by using the Cai-Kim temperature and taking into account an 
infinitesimal time, we get the same result for the horizon entropy as the result obtained in Eq. dSO]). 


III. THE SECOND AND GENERALIZED SECOND LAWS OF THERMODYNAMICS 


On one hand, since cosmos is enclosed by the apparent horizon, it forms a closed system and therefore, the entropy 
of its horizon should increase during the universe expansion meaning that (t^ 


dSA 

dt 


> 0 . 


(56) 


It is called the second law of thermodynamics. Whereas, on the other hand, the total entropy of the closed systems 
should be increased. Since cosmos includes spacetime and its contents, which includes the fluids supporting the 
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geometry of background spacetime, its total entropy consists of two parts including the horizon (Sa) and the conhned 
fluids components (Sin) [80,[8ll- In fact, the generalized second law of thermodynamics states that the rate of the 
total entropy of cosmos including the horizon and confined fluids entropies cannot be negative or briefly [13, [Hj 


(^Sa dSin ^ 
dt dt ~ 


(57) 


Indeed, the total entropy of gravitational systems should meet (l57l) [13, EH. But, here we point to the required 
conditions for satisfying both of the above criterions. 


A. Non-Interacting case 


For the non-interacting cases and while po meets (12511 . by taking a time derivative of the Friedmann equation dHD 
and using the energy-momentum conservation law (Ei) to get the Raychaudhuri equation 

H = + (58) 

Since during the cosmos life < 0 [H, we get 1 — /3 — > 0 leading to H < (I^) while p+p > 0, and 

1 — /3 — < 0 which yields H > for p + p <G. Using Eqs. (1501) and (IMl) to obtain 


rj. dSA ^ At:{p + p) H 

dt 


(59) 


It seems that horizons may satisfy the second law of thermodynamics meaning that the dSA > 0 condition should be 
valid [13, El) 0, [ill ■ lo order to check the validity of the second law of thermodynamics we insert T = — ^(1 + 
into this equation, and get 


dSA _ Sir'^ip + p) 

~dr ~ IP 


(60) 


meaning that the second law of thermodynamics is available for the apparent horizon whiles p+p > 0. This conditions 
leads to lu; > — 1 for the state parameter oj. Moreover, by combining Eqs. (ESI) and EZD with together, we get 


dSA 27ri7 Pd+Pd. 

^ =--jp (1 + 4,— 


(61) 


which means that the second law of thermodynamics is satisfied if 1 -I- > 0. Finally, the second law of 

thermodynamics > 0) is met by the horizon component when po + pd > P + P > 0 are satisfied, 

simultaneously. It is useful to mention here that one can get 

H =-‘i'n{p + p +Pd+Pd), (62) 


by equating Eqs. ED) and ED, which is nothing but the Raychaudhuri equation obtainable by taking time derivative 
from Eq. ED and using (l24ll . Therefore, when p -l-p > 0 and ^ > —{pn + Pd) are available, then p+p + ^ > 
~{pd +Pd) is obtainable which is in agreement with the Raychaudhuri equation (|62|) . For the fluids confined to the 
apparent horizon with total density p, the Gibbs law implies [9^ . 


Ti, 


dSin 

''~dr 


dE dV ^^dp , AttH 

■P— = V^-{p + p)- 


dt 


dt 


Hi ’ 


(63) 


where > 0 is the temperature corresponding to the confined fluids. Now, using (|55)) and V = to get 


dS^n ^ 4Tr{p + p) 

dt £72 ( + H^P 


(64) 


telling us that, for p + p > 0, > 0 is obtainable when 1 -I- -^ < 0 which leads to H < The latter means 

that for the perfect fluids with state parameter ui which either meets the ui < —lor—i < uj conditions, > Q. 
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Additionally, for a prefect fluid with state parameter —1 < w < — the p + p>0 condition is satisfied but < 0. 
Finally, for a prefect fluid with state parameter w which satisfies the — | < w condition the generalized second law 
of thermodynamics > 0) will be satisfied if the p + p > 0 and po + Pd > — H conditions are valid. It 

is useful to mention here that w = — 1 leads to = 0 and = 0 meaning that the generalized second law of 
thermodynamics is marginally satisfied. Moreover, for a more general manner in which oj is not a constant, by using 
the Raychaudhuri equation, we get 


1 + ^ = 1 -47r(p + p) 


1 


H'^{1 - 13) - 


(65) 


Thus, 1 + -^ < 0 leads to 

H^{1 — (3) — < Att{p + p), (66) 

which indicates that > 0. Therefore, if this condition is valid, then the generalized second law of thermodynamics 
will be satisfied. 

For the p + p < 0 case, it is obvious that, from Eq. dsni), ^ < 0. In addition, when H meets the 1 + -^ < 0 
condition, < 0 and thus < 0 meaning that the generalized second law is not satisfied. Briefly, 

for a prefect fluid with ui < —1, the generalized second law is not satisfied. If the Hubble parameter satisfies the 
1 + -^ > 0 condition Eq. del leads to > 0 and therefore, it is legally possible to meet the generalized second 
law of thermodynamics. Using Eq. (l5^ to get 


1 + ^ = 1 -47r(p + p) 


FI2(1 -/3) - 


(67) 


Therefore, the 1 + -^ > 0 condition leads to 

47r(p+ p) < i7^(I —/3) — (68) 

Einally, we can say that if this condition is valid, then > 0 which may lead to satisfy the generalized second law 
of thermodynamics. 

For the horizon entropy of the flat FRW universe supported by a DE candidate with unknown density profile pn, 
we can use Eqs. (HU) and (na) to obtain 

dSA 8n^ip + p) 

sr = 

whenever, it is easy to check that Eqs. (1641) and (IFTl) are also valid in this manner. Bearing Eq. (15^ in mind, H < 0 
leads to p + p > —po — Pd- Similarities with the previous case, in which pD meets Eq. (Ha, are obvious. In fact, 
in order to achieve a more detailed resolution about the validity of generalized second law of thermodynamics, we 
need to know the dependence of either pD or p to the Hubble parameter. We should note again that the horizon 
component satisfies the second law of thermodynamics by (^^ > 0 ) if the pd + pd > and p + p>Q conditions 
are met, which is in agreement with the Raychaudhuri equation (15^ . Moreover, since > 0 is valid when — | < w, 

the generalized second law of thermodynamics > 0 will be available if the — ^ < w and pD + pd > —^ 

conditions are met simultaneously. More studies on the availability of the second law of thermodynamics and its 
generalization needs to know the exact form of p. 


B. Interacting Case 


Eor this case, by using we get again 

dSA _ 27rij 47rpji(l+a;ji) 

dt ^ H ^ ’ 

On one hand, when pd{1 + ujd) > since observationally < 0 [l|, it seems that > 0 is valid everywhere. 
On the other hand, by combining Eqs. dsni), dn and (l23l) , once again we get 

dSA _ 8Tr^{p + p) 


dt 


(71) 
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meaning that > 0 is valid everywhere, ifpD(l+a;D) > —and p+p > 0 are satisfied simultaneously. Therefore, 
the quality of validity of the second law of thermodynamics is similar with the non-interacting case. In addition, 
Eq. (1551) leads to 


+ §^]- VH[p'^ + +Pd)], (72) 

where prime denotes derivative with respect to the Hubble parameter, again. Here, we focus on the model proposed 
by Lime st al. [s^. In this model, wd = — 1 while the vacuum density meets Eq. (I25|) . Substituting into the above 
equation to get 


dS,n Att^p + p) H ^ 3I3H + i 

=-- 4 ^- 

Since H < 0, the second term of RHS of this equation {—VHp'jj) is positive everywhere and therefore, the validity 
of > 0 and thus the generalized second law of thermodynamics depends on the value of the first term of RHS 
(_47 r^+p) jg useful to mention here that for a perfect fluid either obeying w < —1 or —i < w, the Hubble 

parameter meets the H < \ condition leading to 1 -|- -^ < 0 and thus > 0. Moreover, from Eqs. dzni) and (ED 
it is apparent that > 0 when ujd = — I and — I < w, respectively. Therefore, for the flat FRW universe embraced 
a prefect fluid which satisfies — ^ < w, the generalized second law of thermodynamics is convinced. As again, more 
studies on the availability of the second law of thermodynamics and its generalization needs to know the exact form 
of p. 


IV. SUMMARY AND CONCLUDING REMARKS 

Throughout this paper, we considered the FRW universe filled by a DE candidate together a fluid, which is the agent 
of the other possible sources, which may include the baryonic and non-baryonic matters, enclosed by the apparent 
horizon of the flat FRW universe. In continue, we proposed a profile density for the DE candidate which covers 
proposals including CC and dynamic models of DE such as ghost dark energy model, its generalization, the Lima’s 
model and etc. Moreover, by taking into account the Hayward-Kodama definition of the temperature definition of 
apparent horizon as well as the Friedmann equation, we could find the horizon entropy for models in which the DE 
candidate does not interact with the other parts of the cosmos. Our study shows that the DE candidate may modify 
the horizon entropy. We have shown that our formula for entropy (13511 is compatible with previous results about the 
ghost dark energy and its generalization [t^, . Indeed, similar result with (HTl) is reported by authors in ref. [t^ . 

But, our derivation is completely different. Here, we have considered the apparent horizon as the causal bound of the 
system, whiles authors in (76| used the trapping horizon as the causal bound to get the associated horizon entropy. In 
addition, we have generalized our formulation to models in which the DE candidate is an arbitrary unknown function, 
and showed that the DE candidate may modify the horizon entropy (H511 independent of the other parts of cosmos. We 
have also used the Cai-Kim temperature to get the horizon entropy, and found out that the same result for the horizon 
entropy is obtainable if one considers an infinitesimal time in which dV = 0. Thereinafter, we focused on the models 
in which the DE candidate interacts with the other parts of cosmos. We found that the mutual interaction between 
the cosmos sectors may also modify the apparent horizon entropy (1501) . Our study shows that for models in which 
ujd = —1, such as the model proposed by Lima et al. [s^, the mutual interaction between the cosmos sectors does not 
disturb the Bekenstein limit of the horizon entropy. It means that there is no modification to the horizon entropy for 
interacting models with uid = —1 and therefore, Sa = is available in these models. The same as the non-interacting 
case, we tried to get a relation for the horizon entropy in the interacting models by using the Cai-Kim temperature. 
Our study shows that the same result as that of obtained by considering the Hayward-Kodama temperature is available 
for the horizon entropy. Additionally, we pointed to the some required conditions for availability of the second law of 
thermodynamics and its generalization in the interacting and non-interacting models. Our studies show that for the 
non-interacting case, whiles po +Pd > the second law of thermodynamics and its generalization are inevitably 
valid if the state parameter of other parts of the cosmos satisfies the — ^ < w condition. It is because > 0 and 

> 0 are separately valid in this situation. Finally, our study shows that for the interacting case with lod = 

> 0 and > 0 will be met if — i < w and therefore, the generalized second law of thermodynamics will be 
available in an unavoidable way. 
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